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We study a lattice field theory described by two flavors of massless staggered fermions inter- 
acting with U{1) gauge fields in the strong coupling limit. We show that the lattice model has 
a SU (2) X SU (2) X U{1) chiral symmetry and can be used to model the two-flavor QCD chiral 
phase transition in the absence of the anomaly. It is also possible to add a coupling to this model 
which breaks the chiral symmetry to SU{2) x SU (2) and thus mimics the effects of the anomaly 
in two-flavor QCD. We construct an efficient directed loop algorithm to study such a model. We 
show that the chiral phase transition in our model is first order in the absence of the anomaly, 
while it becomes second order with (9(4) exponents when the anomaly is turned on. 
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1. Introduction 

Chiral symmetry plays an important role in determining the low energy physics of QCD. It 
is also expected to play a vital role in determining the behavior of the chiral phase transition that 
occurs at finite temperatures. Renormalization group arguments based on £-expansion techniques 
show that the number of light quark flavors and the strength of the anomaly in QCD can play 
a significant role in determining the order of the transition [[I], |2|]. While a second order phase 
transition is possible for two massless quarks, the presence of more light quarks can introduce 
fluctuations that will force the transition to become first order. If the strength of the anomaly 
is sufficiently weak at the transition, the two-flavor transition could turn into a first order one. 
Although £-expansion techniques are not completely justifiable for the problem at hand, there are 
many examples in which they are known to work sufficiently well. On the other hand there are also 
examples where these techniques fail. 

The effect of the number of light fermion flavors on the chiral phase transition has been a 
subject of extensive research in the past two decades [^, ^ ^, ^ 0]. For recent reviews see [||, ^ p^ ]. 
Unfortunately, the effect of the anomaly on the chiral phase transition has not been studied to a 
similar extent in QCD. Of course it is not possible to tune the anomaly easily in QCD. The subject 



has been analyzed in the context of mean field theory [|11|, |12p and the results are in agreement with 
the £-expansion as expected. Thus, if the strength of anomaly at the phase transition was small 
enough it could be another reason to expect a first order transition in QCD. Lattice simulations 
have shown that the strength of the anomaly may indeed be small at the chiral phase transition 



[13, 14] 



Recently, a new analysis based on a re-summation technique has emerged, which seems to 
show that the results of the £-expansion in two-flavor QCD in the absence of the anomaly may 
not be correct [^]. While the £-expansion forbids a second order transition, the new analysis 
allows it. It would be interesting to find the predicted second order phase transition and its critical 
behavior through Monte Carlo methods. Since two-flavor QCD in the absence of the anomaly has 
an SU (2) x SU (2) x ?7(1) symmetry, the new analysis implies that there is a second order finite 
temperature critical behavior in the appropriate 0(4) x 0(2) sigma model. As far as we know, there 
are no Monte-Carlo studies in the sigma model context which try to address this question from first 
principles. 

In this work we use strongly coupled two flavor lattice QED with massless staggered fermions 
to model the physics of the two-flavor QCD chiral phase transition. Although our model is not 
QCD it has the right chiral symmetries of two-flavor massless QCD and can be studied efficiently 
with cluster algorithms. Our model may be viewed as an alternative to the usual sigma model 
approach where symmetries play an important role. Just like in the sigma model we can add a four 
fermion term to our action which mimics the physics of the anomaly. This allows us to study the 
phase diagram as a function of the temperature and the strength of the anomaly. The similarity 
with the sigma model is only true in an indirect sense. More directly, our model is much closer to 
QCD in spirit since the model is constructed with quark degrees of freedom. Of course the quarks 
are confined into pions and we can rewrite our model as a statistical mechanics of pion world lines. 
Thus, our model offers a new method to explain the qualitative features of pion physics of QCD in 
a simple setting. 
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Figure 1: An illustration of a DPI configuration in 1 + 1 dimensions. 



2. Model 



The action of strongly coupled lattice QED that we study is given by 
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(2.1) 



where x is a lattice site on a four dimensional hyper-cubic lattice with dimensions x L^. The two 
component Grassmann fields and describe the two flavors of quarks and <^aj: represents the 
U{\) gauge field associated to the bond connecting the sites x and x + a (a = 0, 1,2,3 represents the 
temporal and spatial directions respectively. Further, r\a,x are phase factors satisfying the relations 
{^t,xY' = T and (t],-,jc)^ = 1 for / = 1,2,3. We will call the parameter T as the temperature. Thus, 
we have two tunable parameters T and C in our model. 

It is possible to integrate out the gauge fields completely and rewrite the partition function 
as a statistical mechanics of Dimer-Pion-loop-Instanton (DPI) configurations. Each configuration 
is made up of four types of objects: u-dimers 7r"„ =0,1, d-dimers 7r^'„ =0,1, oriented pion- 
dimers = 0, 1, — 1, and instantons 4 = 0,2. A dimer is an object connecting neighboring sites 
and is associated with the corresponding bond (x, a). There are constraints in the configurations 
that follow from the fermionic nature of the microscopic action. These can be mathematically 
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represented as follows: 

h + Y^Km + <« + I <« I + K-a,a + 4-a,a + I <-aA\ = (2.2a) 
a 

iKa - <« + K-a.a " 4-a.a] = (2.2b) 
a 

Y}'^la + T^La,a] = ^ (2-2c) 
a 

The Boltzmann weight of each DPI configuration turns out to be T"'P' where fit is the total number 
of temporal dimers (including all the three dimers) and nj is the total number of instantons. An 
illustration of a DPI configuration is shown in figure [T] 

It is easy to verify that when C = the action is invariant under a ?7 (2) x ?7 (2) chiral symmetry: 

^>,^^U¥,^, »F,„^/?»F,,„ W,^^W,R\ ^,„^^,„L1- (2.3) 

where L,R (2) and Xg and Xg refers to even and odd lattice sites. Since the baryon number is 
gauged, in the strong coupling limit the baryon number is confined. Thus, the relevant symmetry is 
SU{2) X SU{2) X ?7^(1) just like the two flavor QCD. When C / 0, the ?74(1) symmetry is broken, 
which suggests that it induces the effects of the anomaly. For this reason we call the objects that 
arise due to C in DPI configurations, as instantons. Our motivation in this work is to study the 
T — C phase diagram of our model. At low temperature we expect the chiral symmetry to be 
spontaneously broken, while at high temperatures the symmetry is restored. Our intent is to study 
the chiral phase transition as a function of T . We fix Lf = 4 in this work. 



3. Observables 

In the low temperature phase we expect SU (2) x SU {2) xUa{^) symmetry to be broken to 
SU (2) symmetry. In other words 0(4) x 0(2) symmetry is broken to 0(3). In order to study the 
nature of the chiral phase transition we focus on the winding number susceptibilities defined using 
one of the 0(4) currents and the 0(2) currents. 

= I3 I ^ (3.1a) 

= I3 I ^(l<« + <a + l<al) ^ (3-lb) 

We will define = Yimi^ooXv and = limz,^oo Ya- Both F^ and F^ are the squares of the effective 
three dimensional non-singlet and singlet pion decay constants and should be non-zero when the 
corresponding symmetries are broken. On the other hand in the high temperature phase when the 
symmetries are restored we expect the decay constants to vanish. If the chiral phase transition is 
second order we expect 

LY=f{tL'''') (3.2) 

where /(x) is a universal function analytic at x = and v is a critical exponent, where t = T — Tc/Tc 
is the reduced temperature. 



4 



Ejfects of the anomaly on the QCD chiral phase transition 



Shailesh Chandrasekharan 




Figure 3: Plot of F"^ and as a function of T. The data fits well to a power law except very close to Tc. 
The inset shows a two state signal in Yy (as a function of simulation time) at 7" = 2.4675 and L = 48. 



4. Results 

When C = the anomaly is absent. As discussed in the introduction, early studies based on 
the £ -expansion had indicated that the transition must be first order On the other hand recent 
studies have suggested the possibihty of a second order transition with new critical exponents 



[]15[]. So what does our model show? We have studied this question at L, = 4 for lattice sizes 
L = 24,32,48,64 and at different temperatures. In figure ^ we plot LYa (left plot) and LYy (right 
plot) as a function of T for various different values of L. The existence of a second order transition 



can be confirmed if all the curves intersect at a critical temperature as predicted by eq.(3.2). 
Clearly, since the curves do not intersect at a point we conclude that in our model the transition is 
first order. Our data suggests that T^. ~ 2.467(2). 



We have also computed and F^ from the large L extrapolations of Yy and Ya and these 
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Figure 4: Plot of LYv as a function of T for different lattice sizes at C = 0.3. The existence of a second 
order transition is clearly seen since all the curves intersect at Tc = 2.83575(5). The inset shows LYy versus 
(r - T,)L^I'' with V = 0.745 expected for a M 0(4) model. 

results are plotted in figure || as a function of T . Away from tiie critical point both and 
fit well to a power-law of the form A{Tc — Ty . However, the fit gives Tc = 2.4693(3) which is 
too high, confirming again that the transition is indeed first order. In the inset of figure |3| we plot 
evidence for a two state signal at T = 2.4675, a value which is close to the critical temperature. 

Next we switch on the anomaly and study our model at C = 0.3. In this case the symmetry of 
our model reduces to SU (2) x SU (2) which is broken to the diagonal SU (2) at low temperatures. 
Now the finite temperature phase transition can indeed be second order in the universality class of 
the three dimensional 0(4) spin model. This universality class was studied recently in [ p^ ] and it 
was found that v = 0.745(2). In figure ^ we plot YyL as a function of T for values of L exactly like 
figure ^. However, now the curves for different L's do intersect very nicely at a point. In fact we 
can fit our data, close to the critical point to the form YyL = fo + f^tL^/^ + fjtL^^^ , with v = 0.745. 
Such a fit reveals that Tc = 2.83575(10). In the inset of figure § we plot LYy as a function of tL^ /v. 
We see that all of the data collapses nicely to a single function as expected. Thus, we confirm that 
our model shows 0(4) critical behavior at a finite non-zero value of C. 

In summary, thermodynamics of strongly coupled QED with two flavors of staggered fermions 
including a four-fermion coupling C, which models the chiral phase transition in two flavor QCD 
with a tunable anomaly, shows a first order transition at C = and an 0(4) critical behavior at 
C = 0.3. In particular we do not see the second order transition predicted in [jlSl]. We are currently 
trying to locate and study the properties of the tri-critical point present in the C — T plane. 

Part of this work was done in collaboration with D.J.Cecile. For a theoretical discussion of the 
properties of our model and the algorithm used we refer to his contribution to these proceedings. 
This work was supported in part by a DOE grant DE-FG02-05ER41368. 
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